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In this paper we present results of large-scale correlated 
calculations of triplet photoinduced absorption (PA) spec- 
trum of oligomers of poly-(para)phenylenevinylene (PPV) 
containing up to five phenyl rings. In particular, the high- 
energy features in the triplet PA spectrum of oligo-PPVs are 
the focus of this study, which, so far, have not been investi- 
gated theoretically, or experimentally. The calculations were 
performed using the Pariser-Parr-Pople (PPP) model Hamil- 
tonian, and many-body effects were taken into account by 
means of multi-reference singles-doubles configuration inter- 
action procedure (MRSDCI), without neglecting any molecu- 
lar orbitals. The computed triplet PA spectrum of oligo-PPVs 
exhibits rich structure consisting of alternating peaks of high 
and low intensities. The predicted higher energy features of 
the triplet spectrum can be tested in future experiments. Ad- 
ditionally, theoretical estimates of exciton binding energy are 
also presented. 
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I. INTRODUCTION 

Photoluminescent conjugated poly- 

mer poly- (para)phenyleneviny lane (PPV) is one of the 
most promising candidates for the optoelectronic devices, 
of the next generation, such as light-emitting diodes,El 
and lasers.EI Because of their quasi one-dimensional char- 
acter, the electrons(e) and holes(h) in conjugated poly- 
mers have a high recombination probability, -leading to 
their large photoluminescence (PL) efhciency.El However, 
in electroluminescence experiments, a significant fraction 
of the e-h recombination events lead to ths,formation of 
triplet excitons, which are nonluminescentE In addition, 
it is believed that through various other processes such 
as intersystem crossings and fission,B the singlet exci- 
tons in these materials transform into the triplet ones. 
Therefore, in order to fully characterize the PL prop- 
erties of conjugated polymers such as PPV, a thorough 
understanding of their triplet excited states is required. 
A good understanding of the triplet states is fruitful from 
another vantage point as well. The energy difference be- 
tween the triplet and singlet states defines the exchange 
energy, and consequently, quantifies the strength of e-e 
correlations in the system. Therefore, a deeper under- 
standing of the triplet states will also enhance our un- 
derstanding of electron correlation effects in conjugated 
poWmers. Recently Oster|backa et al.p Romanovskii et 
al.Jj and Monkman et al.El have measured the energies 
of the first singlet (^i) and triplet (Ti) excited states of 
various conjugated polymers. Triplet excited states with 



energies higher than Ti are frequently probed through 
the photoinduced -absorption (PA) spectroscopy in the 
triplet manifold.&Q In these experiments first the Ti is 
populated, and then the optical transitions are induced 
to higher-energy triplet states. Assuming a planar config- 
uration for oligo-PPVs corresponding to symmetry group 
C2/1, their one- and two-photon states will belong, respec- 
tively, to the Bu, and the Ag irreducible representations. 
The ground state of these oligomers is l^Ag, while Si cor- 
responds to l^Bu- In the triplet manifold, Ti corresponds 
to the l^Bu state, while the first state with a strong 
dipole coupling to l^B^, is called rrt'Ag. Therefore, in 
the triplet PA spectrum of oligo-PPVs, m^Ag would show 
up as the first peak, while any higher ^v4g-type states 
with significant dipole coupling to l^B^, would lead to 
high energy features in the spectrum. Therefore, through 
triplet PA spectroscopy, one can probe the triplet two- 
photon states (^An) of ohgo-PPVs. jHowever, so far all 
the theoretical,ErU and experimentaJj^El investigations of 
triplet PA in PPVs have focused on the first rrr'Ag state, 
providing no information on the higher-energy "^A^-type 
states. Considering recent PA measurements of high- 
energy two-photon states in the singlet manifold {^Ag) 
of the PPVs which have yielded a lot of exciting infor- 
mation on the nature of those statesEl, we believe that a 
similar analysis should be extended to the high-energy 
■^^g-type states of this material. Therefore, to fill that 
that void, we decided to undertake a systematic large- 
scale correlated calculation of the triplet PA spectrum of 
oligo-PPVs which will target the high-energy features at 
the same level of computational accuracy as the m^Ag 
peak. However, compared to the singlet- manifold re- 
sponsible for the singlet PA , the correlation effects are 
stronger in the triplet- Ag manifold because here, in ad- 
dition to the orbital excitations, spin excitations are also 
involved. Therefore, it is mandatory that the correlated 
calculations employ large basis sets in order to describe 
the high-energy '^Ag peaks accurately. Since these are the 
first calculations of the kind, for the purpose of bench- 
marking we also performed calculations on a number of 
other excited states of oligo-PPVs, and compare them 
to the works of other authors. As a matter of fact our 
calculations on the triplet PA spectrum indicate the exis- 
tence of three peaks beyond m^Ag exhibiting rich struc- 
ture with alternating high and low intensities. We hope 
that this work will encourage experimental measurements 
of the high-energy region in the triplet PA spectrum of 
PPVs, just like similar recent measurements of their sin- 
glet PA spectrum.Ll 

The remainder of this paper is organized as follows. In 
section || we briefly describe the correlation approaches 
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used to perform the calculations in the present work. 
Next in section [II we present and discuss the calculated 
trip let PA spectrum of oligo-PPVs. Finally, in section 
|lV| we summarize our conclusions, and discuss possible 
directions for future work. 



II. METHODOLOGY 

Henceforth any oligomer of PPV containing n phenyl 
rings will be referred to as PPV-n, making trans-stiVoene 
synonymous with PPV-2. We performed correlated cal- 
culations on the oligomers PPV-2 to PPV-5, assuming 
completely planar geometries and the symmetry group 
C2/i. rfirom among a variety of model Hamiltonians avail- 
able,t3 we chose the Pariser-Parr-Pople (PPP) model for 
the present calculations, which describes the 7r-electron 
dynamics of conjugated systems in terms of a minimal 
basis set, and a small number of parameters. The PPP 
Hamiltonian reads 



(1) 



i<j 



where (ij) implies nearest neighbors, cj^ creates an elec- 
tron of spin (T on the Pz orbital of carbon atom i, 
Uia- — c\^Cicr is the number of electrons with spin cr, 
and rii — riia is the total number of electrons on 
atom i. The parameters U and Vij are the on-site and 
long-range Coulomb interactions, respectively, while tij 
is the nearest-neighbor one-electron hopping matrix ele- 
ment. The Coulontbiinteractions are parametrized as per 
Ohno relationship,ll3 
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where Rij is the distance (in A) between the sites i and 
j, and is the dielectric constant to account for screen- 
ing. Since the results obtained with the PPP Hamilto- 
nian depend substantially on the choice of the Coulomb 
parameters, we tried two sets: (a) "standard parameters" 
with U = 11.13 eV and = 1.0, and (b) "screened pa- 



rameters" with U = 8.0 eV and Ki 



1.0, and Ki 



2.0 



with i ^ j. The screened parameters., originally in- 
troduced by Chandross and Mazumdam in their study 
of optical absorption in PPV, are known to give better 
agreement with experiments on excitation energies, pos- 
sibly by taking interchain screening effects into account. 
For the nearest-neighbor hopping matrix elements tij , we 
chose —2.4 eV for the phenylene ring, and —2.2 eV and 
—2.6 eV, respectively, for the single and double bonds of 
the vinylene linkage. As far as the bond lengths are con- 
cerned, in the phenyl rings they were taken to be 1.4 A. 
In the vinylene linkage the single (double) bond lengths 
were taken to be 1.54 A (1.33 1). 



Starting with the restricted Hartree-Fock (HF), 
depending on the size of the oligomer, the elec- 
tron correlation effects were included either by the 
quadruple configuration-interaction (QCI) method or by 
multi-reference singles-doubles configuration-interaction 
(MRSDCI) approach. Botb_the approaches are well doc- 
umented in the literature. Ej The accuracy of the MRS- 
DCI calculations depends on the choice of the reference 
configurations, and their number [Nref)- The reference 
configurations are chosen according to their contribution 
to the many-particle wave function of the targeted state, 
while the convergence of our results with respect to iVre/ 
will be examined in the next section. To reduce the size 
of the Hamiltonian matrices, we made full use of the spin 
and the C2h point-group symmetry. The triplet PA spec- 
trum was calculated using the oscillator strengths com- 
puted under the electric dipole approximation. The finite 
lifetime of the states was taken into account by means of 
a linewidth parameter F. Further technical details can be 
found in our recent work, where we used this symmetry- 
adapted CI methodology to study the low-hting excited 
states in phenyl-substituted polyacetylenes.tj 



III. RESULTS AND DISCUSSION 

The correlation method, N^ef, and the total number 
of configurations {N total) involved in the CI calculations 
of excited states of various oligo-PPVs are detailed in 
table |. From the sizes of the CI matrices diagonalized 
it is clear that these calculations were very large scale, 
and therefore, they must include the infiuence of electron 
correlation effects to a very high accuracy. The conver- 
gence of our MRSDCI results with respect to Nref^ is 
demonstrated in Fig. |^. It is clear from the results that 
good convergence in the excited state energies has been 
achieved in all the cases by the time fifteen most impor- 
tant configurations {N^ef = 15) have been included in the 
reference space in various MRSDCI calculations, leading 
us to believe that the excitation energies presented here 
are accurate to within a few hundredths of an eV, for the 
given values of PPP parameters. 

The noteworthy aspect of the preseaL-t^culations as 
compared to those of other authorgl3'E3Eil is that no 
MOs have been discarded, from the correlated calcula- 
tions. Lavrentiev et al.JHI with the purpose of obtain- 
ing results for infinite PPV, performed density-matrix 
renormalization group calculations in conjunction with 
a PPP-Hamiltonian-based two-state model consisting of 
just the highest-occupied and the lowest unoccupied MOs 
(HOMO and LUMO) of each phenylene and vinylene 
unit. Thus, e.g., for PPV-5 which has 38 MOs in to- 
tal, the. number of active MOs considered by Lavrentipz. 
et al.EJ was just 18. Additionally, Lavrentiev et al.EJ 
used a different set of PPP jparameters as compared to 
us. Similarly Gartstein et al.E3 simplified the PPP model 
to obtain analytical results for the energy of l^B^ state. 
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but only for infinite PPV. Rohlfing et al.t2l performed 
first-principles calculations on both the sipglet and triplet 
states of infinite PPV. Beljonne et al.JlJ performed cal- 
culations on the same oligomers of PPV as considered by 
us, utilizing an intermediate neglect of differential overlap 
(INDO) model Hamiltonian, also within the framework 
of the MRSDCI procedure. However, they used only four 
reference configurations {Nref = 4), and considered sin- 
gle and double excitations from a set of five highest oc- 
cupied-JT orbitals into the five lowest unoccupied tt* or- 
bitals.t2l Thus, N^ef as well as the number of active MOs 
used by Beljonne et al.tj is much smaller compared to 
our calculations. Therefore, to benchmark our compu- 
tational approach and the choice of PPP Coulomb pa- 
rameters, we compare our MRSDCI results on various 
low-lying excited states of PPV-n [n = 2 — 5) to those of 
other authors and experiments, in table ||. 

Inspection of table |l| reveals the following trends: (i) 
Our MRSDCI excitation energies computed with the 
screened Coulomb parameters are lower, and in bet- 
ter agreement with the experiments, than those com- 
puted with the standard parameters, (ii) For E(\^By) 
our calculations predict complete saturation, in excellent 
qualitative agreement with the highly-localized picture 
of the l^Bu state, verified recently in the experiment of 
Osterbacka et al.p reporting that \^Bu was a Frenkel ex- 
citon, with a radius of just 3.2 A. However, neither our 
calculations, nor those of other authors are close to the 
experimental value of E{l^Bu) ~ 1.5 eV for the infinite 
PPV.qQ (iii) As far as comparison with other authors 
is concerned, our screened parameter values of E{l^Bu) 
are in similar asceement with the experiments as those of 
Beljonne et al.JlJ while thi2v are somewhat higher than 
those of Lavrentiev et al.Eil whose E{l^Bu) for PPV- 
5 is lower than the experimental value. However, our 
screened parameter values of E{m^Ag) — E{l^Bu) are 
consistently in much better agreement witii experiments 
as compared to those of Beljonne et al.ll3 and Lavren- 
tiev et al.El This is very encouraging because, in order to 
compute the triplet PA spectrum accurately, one needs 
to have an accurate representation of states in the '^Ag 
manifold, starting with the m^Ag state. 

To further benchmark our approach, we use E{m^Ag) 
computed above to obtain a theoretical estimate of the 
binding energy of the exciton. Shimoi and Mazum- 
dar have argued that, because m^Ag state is also an 
exciton, E{m^Ag) — E{l^Bu) is a good lower-limit es- 
timate, £'f,(min), for the exciton binding energyHj From 
our calculations, for PPV-5 one gets i?f,(min) — 0.68(1.7) 
eV with screened (standard) parameters. Keeping in 
mind that generally our screened parameter based cal- 
culations were in good agreement with the experiments, 
we take the value i?b(min) = 0.68 eV to be the correct 
theoretical estimate. Thus our theoretical estimate of 
£'fc(min) = 0.68 eV obtained from PPV-5 is in very good 
agreement with the experimental value of -Bb(min) = 0.55 
eV estimated by Osterbacka et al.ul for infinite PPV. 



Next we turn our attention to the main focus of this 
work — the triplet PA spectrum of oligo-PPVs. The cal- 
culated spectrum of PPV-4 is presented in Fig. ||, and 
of PPV-5 in Fig. jsj The spectra have been computed 
both with the standard parameters (Figs. ||a and |^), 
as well as with the screened parameters (Figs. ||b and 
. Thus the four presented spectra will help us under- 
stand the influence of the following two factors on the 
triplet PA: (i) the conjugation length (ii) the strength of 
Coulomb interactions. 

The following trends emerge from the spectra pre- 
sented in Figs. H and |[ (i) Irrespective of the Coulomb 
parameters used, both for PPV-4 and PPV-5 there are 
four major structures (labeled I, II, III, and IV) of alter- 
nating high and low intensities, (ii) Irrespective of the 
Coulomb parameters, the relative intensities of the high 
energy peaks (II, III, and IV), as compared to peak I, are 
increasing with increasing conjugation length (iii) For a 
given oligomer, the relative intensities of the high-energy 
peaks, compared to peak I, are higher when computed 
using screened parameters, as compared to the ones 
calculated using the standard parameters. Although, 
standard-parameter-based calculations yield rather small 
relative intensities for peaks II of PPV-4 and PPV-5, 
however, the same parameters lead to significant intensi- 
ties of other high-energy peaks (III and IV). Thus, based 
upon these results, we arrive at the following conclusions 
regarding triplet PA in oligo-PPVs: (i) The intensities of 
peak III and beyond are strong enough to be detected ex- 
perimentally in rather small oligomers such as PPV-4 and 
PPV-5. (ii) The intensity of peak II for small oligomers 
appears to be strongly dependent on the strength of 
Coulomb interactions. This peak will be detectable for 
these oligomers only if screened parameters, perhaps due 
to interchain screening, describe the strength of Coulomb 
interactions, (iii) For relatively longer oligomers, even 
peak II will acquire significant intensity both for standard 
and screened parameters, as exemplified by the rapid in- 
crease in its intensity with size. 

Having discussed the general qualitative features of 
the computed triplet PA spectrum, we next examine 
the nature of many-particle ^Ag states contributing to 
various peaks. Here we limit our discussion to the 
specific case of the spectrum of PPV-5 computed with 
screened parameters as presented in Fig. |^. The rea- 
sons behind this choice are obvious — PPV-5 is the longest 
oligomer studied, and the values of the energy gaps 
E{m^Ag)—E{l^Bu) with screened parameters for various 
oligomers were in very good agreement with the experi- 
ments. However, we did examine the nature of the peaks 
in the triplet spectra presented in the remaining three fig- 
ures (Figs. |2^, and ^), and all of them were found to 
be qualitatively very similar to the states corresponding 
to the peaks in Fig. Thus we believe that our dis- 
cussion of triplet PA spectrum of PPV-5, computed with 
screened parameters, represents the general features of 
the triplet PA spectra of oligo-PPVs. The first peak in 
Fig. ||b, occurring close to 1.9 eV, is the strongest in the 
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spectrum, and corresponds to the l^Bu m^Ag absorp- 
tion. Compared to the Hartree-Fock configuration, the 
many-particle wave function of m^Ag is a hnear combi- 
nation consisting mainly of singly excited configuration 
(H-1 L) and its charge conjugated (c.c.) counterpart 
(H ^ L+1), where H=HOMO and L=LUMO. Feature 
II located near 2.7 eV has a smaller intensity compared 
to I, and corresponds to two closely-spaced '^Ag states 
also consisting predominantly of single excitations (H — ^ 
L-l-3) -I- c.c, and (H-1 L+2) + c.c. The intensity of the 
next feature (III) is much stronger compared to that of II, 
and is composed of states with main contributions from 
double excitations. The dominant peak of III, located 
near 3.9 eV, corresponds to a ^Ag state composed pre- 
dominantly of double excitation (H-1 ^ L ; H — > L-|-l). 
Feature IV has weaker intensity compared to HI, and 
its main peak occurs close to 5.1 eV. Similar to III, this 
feature also corresponds to states composed of doubly 
excited configurations. The state constituting the main 
peak of IV is predominantly composed of (H — > L+2 ; H-2 
— > L), and other double excitations. In the infinite poly- 
mer limit, we speculate that the single-excitation-based 
features I and II will merge to form the first absorption 
band, while the double-excitation-based features III and 
IV would merge to form the second band. Thus, for 
infinite PPV, we expect the triplet PA spectrum to be 
similar in appearance to its singlet counterpart, which 
also exhibits two major bands PAl and PA2Eil 

Experimental measurement of tuplet PA peaks beyond 
nv^Ag (peak I) is not an easy task.E3 However, if it is pos- 
sible to conquer the associated practical difficulties, we 
believe that the measurements of the high-energy region 
of the triplet PA will shed light on some important ex- 
cited states of PPV, thereby further clarifying the role of 
electron correlation effects in this substance. 



IV. CONCLUSIONS AND FUTURE DIRECTIONS 

In conclusion, very large-scale correlated calculations 
using the PPP Hamiltonian were performed on the oligo- 
PPVs, and their triplet PA spectrum extending well into 
the high-energy region was computed. The computed 
spectrum predicts rich structure beyond m^Ag, consist- 
ing of peaks with alternating high and low intensities. 
The peaks were shown to result from ■^Ag-type states 
composed predominantly of configurations which, besides 
the spin excitation, also exhibit either single or dou- 
ble charge excitations with respect to the ground state. 
Therefore, it will be interesting to probe the higher- 
energy triplet PA spectrum both of oligomers, as well as 
of infinite PPV experimentally, as it will lead to further 
insights into the low-lying excited states of this techno- 
logically important material. Analogous to the singlet 
manifold, it will then also be of interest to explore the 
existence of triplet biexcitons. 

Although the electron-correlation effects were included 



to a large accuracy in the present calculations, yet, be- 
cause of the use of a rigid-lattice model, the influence of 
electron-lattice coupling was^-cpmpletely ignored. How- 
ever, recently Barford et al.ES have demonstrated that 
the contribution of electron-lattice relaxation can be im- 
portant on the low-lying correlated triplet states. They 
demonstrated that for polyenes l^B^ state can relax by 
a few tenths of an eV thereby increasiag the correspond- 
ing E{mPAg) — E{l^Bu) energy gapB3 Therefore, it will 
be of considerable interest to perform a similar calcu- 
lation for oligo-PPVs, and compute the contribution of 
electron-phonon coupling on the triplet excitation ener- 
gies of these materials. These aspects as well as phospho- 
rescence and triplet electroluminescence in conjugated 
polymers will be studied theoretically in future publi- 
cations. 
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FIG. 2. Photoinduced absorption spectrum from V^Bu 
state in PPV-4: (a) computed with the standard parameters 

(b) computed with the screened parameters. A linewidth pa- 
rameter r = 0.1 eV was used. E is the energy of the incident 
photon. 



6 




FIG. 3. Photoinduced absorption spectrum from l^Bu 
state in PPV-5: (a) computed with the standard parameters 

(b) computed with the screened parameters. A Unewidth pa- 
rameter r = 0.1 eV was used. E is the energy of the incident 
photon. 
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TABLE I. For difFerent states of various oligomers, the 
number of reference configurations (Nref) and tiie total num- 
ber of configurations (Ntotai) involved in the MRSDCI (or 
QCI, where indicated) calculations. 



Oligomer l''Ag 



Nref Ntotai Nref Ntotai Nref Ntotai Nref Njot. 



PPV-2 
PPV-3 
PPV-4 
PPV-5 



1" 49495 
r 2003907 
35 1627923 
25 3257930 



r 82423 
1" 3416371 
36 1503239 
25 2537197 



1" 144543" 
65 1212746 
35 2136547 
25 3217908 



72 112193 
65 1321885 
35 2766111 
20 3314575 



" QCI method 



TABLE II. Comparison of our results with the standard 
parameters (This work-1) and screened parameters (This 
work-2) with those of other authors and experiments. All 
the energies are in eV. 



Work 



PPV-2 PPV-3 PPV-4 PPV-5 



1M„ -I'Bu 



iMo -i^Bu 



rB„ -m?A, 



This work-1 
This work-2 ^ 
Beljonne et al.E£ 
Lavrentiev et al, 
Exp. 

This work-1 
This work-2 |- 
Beljonne et al.c£ 
Lavrentiev et al. 
This work-1 
This work-2 |- 
Beljonne et al.t: 
Lavrentiev et al, 
Exp. 
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1.95 
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3.74 
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2.92 


2.90 


3.25 


2.55 


2.14 


1.93 


4.15 


3.26 


2.75 


2.46 


3.63 


3.04 


2.61 


2.34 


3.29'= 


2.30"^ 


1.95^* 


1.80'* 



° Refl 
^ Ref.E 
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